Magnetoexcitons and optical absorption of bilayer-structured topological insulators 
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The optical absorption properties of magnetoexcitons in topological insulator bilayers under a 
strong magnetic field are theoretically studied. A general analytical formula of optical absorption 
selection rule is obtained in the noninteracting as well as Coulomb intra-Landau-level interacting 
cases, which remarkably helps to interpret the resonant peaks in absorption spectroscopy and the 
corresponding formation of Dirac-type magnetoexcitons. We also discuss the optical absorption 
spectroscopy of magnetoexcitons in the presence of inter-Landau-level Coulomb interaction, which 
becomes more complex. We hope our results can be detected in the future magneto-optical experi- 
ments. 
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The bilayer n-p systems, comprising electrons from the 
n layer and holes from the p layer, have been the subject 
of recent theoretical and experimental investigations in 
low-dimensional condensed matter physics. The coupled 
quantum wells [lj-ll] and layered graphene 0-[H| are of 
particular interest in connection with the possibility of 
the Bose-Einstein condensation and superfluidity of indi- 
rect excitons or electron-hole pairs. On the other hand, 
as a new state of quantum matter, topological insulators 
(TIs) p^ - [20j are now being in intensive study. They are 
characterized by a full insulating gap in the bulk and 
topologically protected gapless edge or surface states in 
low dimensions, in which by their single-Dirac-point na- 
ture it is highly desirable to find stable and intriguing 
electron-hole excitations. 

With the rapid progress in advanced nanotechnology, it 
is possible to fabricate a bilayer n-p system with TIs. We 
call this system as "topological insulator bilayer (TIB)", 
which comprises two TI thin films separated by a di- 
electric barrier. The schematic scheme is shown in Fig. 
QJa). Bi2Se3-family materials are one of the best can- 
didates for their large bulk-gap width (21rl23j . In order 
to form steady excitons or magnetoexcitons, the dielec- 
tric spacer can be chosen from those with small dielectric 
constant, such as AI2O3 and SiC>2, on which, remarkably, 
high-quality TI quantum well thin films have been now 
successively grown j24U27| . By doping or applying ex- 
ternal gates, electron and hole carriers will form in both 
topological insulator thin films, which behave like mass- 
less Dirac particles. However, these quasiparticles can 
not form excitons because there is no gap opening. To 
produce a gap we apply a strong perpendicular magnetic 
field on the TIB. In this case, the Dirac-type energy spec- 
trum of these quasiparticles becomes discrete by forming 
Landau levels (LLs), which results in the possible forma- 
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FIG. 1: (Color online) (a) Scheme of TIB consisting of two 
Be2Se3 thin films separated by dielectric spacer. Electron 
(hole) carriers are induced by n(p)-type doping or applied 
external gates. Indirect magnetoexcitons can be formed in 
presence of strong external magnetic field, (b) and (c) are 
schematic representations of the chemical potential location, 
(b) represents the fully occupied case and (c) represents the 
partially occupied case. 



tion of magnetoexcitons in this system. 

The optical absorption spectroscopy analysis is an in- 
structive method in studying the properties of magne- 
toexcitons, because it can provide the knowledge of mag- 
netoexciton's energies and wavefunctions at P=0, where 
P is the magnetoexciton momentum. Many efforts have 
been devoted to exploring the optical absorption proper- 
ties of magnetoexciton in semiconductor quantum wells 
by using the magneto-optical measurements [H, __!_]. 
More recently, the optical characterization of Bi2Ses in 
a magnetic field has been experimentally studied [3Cj . 
Based on this, the main purpose of the present pa- 
per is to study the magnetoexciton's optical absorption 
spectroscopy in TIBs. Note that the magneto-optical 
response of Bi2Se3 thin film grown on dielectric sub- 
strate has recently been experimentally measured [27| . 
Through applying a strong magnetic field, the Coulomb 
interaction (CI) between electrons in the upper TI film 
and holes in the lower TI film can be thought as a pertur- 
bation compared to the energy difference between LLS. 
For convenience, we divide the CI perturbation into the 
intralevel [see the following text and Eq. ((SJ)] and inter- 
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level parts. Firstly, we consider a simple case in which 
only the intralevel CI is introduced. In this case, we de- 
duce a general analytical formula for magnetoexciton's 
optical absorption selection rule. According to this selec- 
tion rule, one can clearly point out which magnetoexciton 
corresponds to the specific resonant peaks in absorption 
spectroscopy. Then, we study the case in which inter- 
level component of CI is also included. In this case, the 
optical absorption spectroscopy becomes more complex 
and intriguing by the mixing of the excitations. 

We start with the effective Hamiltonian of TIB system 
H=H +U(r), where 



Ho = v F a e ■ (z x 7r e ) - v F a h ■ (z x n h ) 



(1) 



is the free electro n-hole part of TIB, while 
U(r)=— e 2 /e-y/|r| 2 + d 2 is the CI between the pair 
of electron and hole with r=r e — r;, the relative coordi- 
nate in the x-y plane and d the spacer thickness. Here, 
r e(h) represents the position of the electron (hole), and 
7r e (fe) =P e (fc)±eA eW /c=-id/<9r e(h )±eA e (/ l )/c denotes 
the in-plane momentum of the electron (hole), where the 
gauge is chosen as A e ( ft )=B(0, x e ^, 0) in the following 
calculations, vf is the Fermi velocity (~3xl0 5 m/s for 
Bi2Se3-family materials [3l|, HH), z is the unit vector 
normal to the surface, and <r e = a ® I (<Jh= I <£> cr) 
describes the spin operator acting on the electron (hole) , 
in which a denotes a vector of Pauli matrices and I is 
the 2x2 identity matrix. 

The eigenstates of Hamiltonian Hq have the form 



ipp (R, r) = cxp 



i[P+— [Bxr 
2c 



R 



* (r-ro) , (2) 



where R=(r e +r/ 1 )/2, r=r e — r/,, and rQ—l B (BxP) with 
Ib being the magnetic length. For an electron in LL n + 
and a hole in LL n_ , the four-component wave functions 
for the relative coordinate are given by Q 



*« + ,n_(r) = |ra + ,n- 



|n + l-l,|n-| 



-l(r) 



where 



=* r. 



.(r)= 



|n + |-l,|n_|(r) 

« sgn(n+) $|„ + |,|„_|_i(r) 

^sgn(n + )-sgn(n_) $|n+|)|n _ |(r) J 



(3) 

\ 



with i z =Hi— n2, n_ = min(ni, 712), and 

<5 (Z z )=sgn(/ 2 ) iz — >1 for ? z =0. The corresponding 
LLs are given by 



V2- 



sgn(n+) y/\n+\ -sgn(n_ ) y/\n- 1] . (4) 



Taking CI as a perturbation in the first order and only 
considering its intralevel component, the energy disper- 
sion of magnetoexcitons can be easily obtained as 



n_j_ ,n_ 



|£7(r+r )|* r . 



(5) 
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FIG. 2: (Color online) Energy dispersion of the first five mag- 
netoexciton LLs for the fully (left panels) and partially (right 
panels) occupied cases. ,8=0.9 and /i=f . The spacer thickness 
d=0.2l B for (a) and (b), while d=2.0l B for (c) and (d). 



However, to take into account the interlevel CI, we 
should perform diagonalization of the full Hamiltonian 
for Coulomb interacting carriers in some basis of magne- 
toexciton states x I'n + ,n_(r). To obtain eigenvalues of the 
Hamiltonian "H, we need to solve the following equation: 



0= det 6 n+in , + 5 n _ >nL (Ei°l n _ -E) 
n'_ \U (r-j|SxP) |*„ + , n _ 



(6) 



The location of the chemical potential will determine the 
possible LL indices for electrons and holes. For conve- 
nience, in this paper we use the notation fi for the highest 
filled LL and only consider two cases: (i) Electron LLs 
with n+>/i are unoccupied and hole LLs with are 
fully occupied and (ii) electron LLs with n + >[i are unoc- 
cupied and hole LLs with n_</i are fully occupied, while 
the LL at [i is partially occupied. These two cases are 
schematically represented in Fig. [IJb) and 1(c), respec- 
tively. To distinguish these two cases, in the following 
discussion we call the first case as fully occupied case 
and the second one as partially occupied case. 

We take /i=l as an example to numerically solve Eq. 
([6]) by employing five electron and five hole levels. In our 
calculations we use the notation /3= {e 2 /ds) /{vf/Ib) = 
e 2 / (cvf) to describe the relative CI strength. Figures 
[2ja) and [2fb) plot the first five energy levels as a func- 
tion of P for fully and partially occupied cases, respec- 
tively. Here, the spacer thickness is chosen to be d=Q.2lB- 
With increasing d, the CI effect turns to become faint. 
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As a result, the magnetoexciton's energy levels become 
smoother. This can be clearly seen from Figs. [5Jc) and 
Eljd), where the spacer thickness d=2Ms is used for com- 
parison with cases of <1=0.21b- 

Let us now turn to discuss the optical absorption prop- 
erties of magnetoexcitons with P = 0. The magneto- 
optical experiments provides a good method to detect the 
magnetoexcitons since the particle-hole excitation energy 
could be determined directly by the resonant peaks in the 
optical conductivity. Up to this stage, the LLs have been 
assumed to be sharp, and thus from Fermi's golden rule, 
the optical absorption for photons of frequency lo yields 
a sum of delta functions in energy, 




li=0 - #-<>•"•-•' • 

MM 



(b) 



JUL 



1~* — 

=0- "•"O^f-*— 
-1 » > o-»- 

MM 



2.4 



0.8 



1.6 



2.4 



FIG. 3: (Color online) The optical absorption curves with 
fi—0 for (a) fully and (b) partially occupied cases. The spacer 
thickness is set to cI=0.21b- The corresponding schematic 
sketch of magnetoexcitons are plotted in the insets. 



(a\-A- v|0) 

c 



5(e^-fuv), (7) 



where the coefficients c; (i=l, ■ ■■ , 4) are defined as 



where a is the set of quantum numbers describing a 
particle-hole excitation and |0) = |/x, /i) is the ground state 
in absence of particle-hole excitations. v—dHo/fid'k is 
the velocity operator for non-interacting electrons (we 
neglect the renormalization of the Fermi velocity due to 
interactions). A is the vector potential. For linear polar- 
ized light A=Ax, one easily deduce v x ccm^ — (a^-af i )/2. 

For the purpose of resolving the spatial functional 
forms and reflecting realistic experimental conditions 
one should either assume that the LLs are broad- 
ened or that the magnetoexciton has a finite lifetime. 
We choose, therefore, a Lorentzian type broadening 



7/2 



with linewidth 7 in following 



( e ^'"^) 2 +7 2 /4 

calculations. Substituting it into Eq. (O, one gets 



ci= 

C2 = 
C3 = 
C 4 = 



^_^\-sgn(n-)_|_^l+sgn( ( u)/_^sgn(n+)-Bgn(n-) 



^sgn(n + )4_,U-sgn(/i) /_^\sgn(n + )-sgn(n_) 



(_j) 1 - s g n (M)4_(_j)l+sgn(n + ) 



c. 



•l+sgn(>i) 



+*(-*) 



-sgn(n_) 



c, 
c, 

(11) 



with C=(V2) Sn +- 0+Sn -- 0+2S » '°~ 4 . Our optical absorption 
selection rule formula (|10|) proves to be very useful when 
analyzing magnetoexciton's absorption resonant peaks in 
spectroscopy of TIB systems. For the fully occupied case 
with fi^O, the selection rule can further be simplified as 



(n+,n-\m x \n,fi) = c 4< 5 n+iAI+ i<5„_ 



(12) 



R(lo)(x 



7 /2 



(8) 



with |(a|m-|0)| = E» + > M Tdr fOT ful1 occu " 

pation and En En_^ J ^ r l° r partial occupation. 
The parameter T is defined as 



T =\ C n+,n- (n+,n- \m x \n,fj) I , 



(9) 



where C" n _ is the projection of magnetoexciton state 
a on the basis state |n + ,n_). 

Let us now consider a simple case in which only the 
intralevel CI is considered. This case is similar to that 
without CI, in which the magnetoexciton state a is one 



special basis state and C" iT( =5 a ,( n+) „ )• So by sub- 
stituting Eq. ([3]) into n_|m~|/x, we can imme- 
diately obtain the optical absorption selection rule for 
magnetoexcitons as 

(n+,n-\m~\fj,,fi) (10) 

=Cl<J|„ + |,| Al |5|„_|,| A( |_i+C25| n+ | > | A( |_ 1 5| n _| i | Al | 

+C3<Vt-l>H V-UmI+i+^V+I.H+iV- l>H ' 



The present general optical absorption selection rule 
proves to be very useful when analyzing magnetoexci- 
ton's absorption resonant peaks in spectroscopy of TIB 
systems. 

For simplicity, let us first consider the fully occupied 
case with /^=0, i.e., electron LLs n + >0 arc unoccupied 
and hole LLs n_^0 are fully occupied. The ground 
state is |0)=|0, 0). According to the selection rule (|12l) . 
the non-zero elements of n_ \m~ |0, 0) only occurs at 
n+=l while n_=0. That means there is only one reso- 
nant peak to appear in the absorption spectrum, which 
corresponds to the energy of magnetoexciton state |1, 0) 
(see Fig. EJa)). For comparison, we also study the par- 
tially occupied case with /i=0. According to the selection 
rule (ITU|) , (n+ , n_ I m~ 1 , 0) ^0 only when n+=l , n_ =0 or 



=0, n_ 



T. Thus in this case there are two resonant 



peaks occurring in the absorption spectrum as shown in 
(see Fig. [3jb)), which correspond to the energies of two 
magnetoexciton states |1,0) and |0, — 1). By increasing 
d (equally decreasing the CI), these two magnetoexciton 
levels 1 1 , 0) and |0, — 1) turn to be degenerate, which re- 
sults in the incorporation of these two resonant peaks to 
be a single one at d—too. Schematic sketch of the corre- 
sponding magnetoexciton state obeying the selection rule 
(IT21 are plotted in the insets of Fig. [3J 
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faint and that at |4, 1) disappears (blue dashed curve s 



FIG. 4: (Color online) Same as Fig. [3jb) except fj,=l. The 
green dotted curve represents the corresponding optical ab- 
sorption spectrum at the limit of d— >oo. 



To ulteriorly see the difference of absorption proper- 
ties between partially and fully occupied cases, let us 
now consider the case of //=1, i.e., the ground state is 
|1, 1). The selection rule for fully occupied case (Eq.(fT2|)) 
with /i— 1 promises that there is only one resonant peak 
occurs, which corresponds to the formation of magne- 
toexciton state |2, 1) by absorbing a photon quanta u>. 
In the limit d— >oo, lu— >v / 2(v / 2— 1)=0.59. For partially 
occupied case, however, there are three non-zero ele- 
ments according to selection rule flTU|) , which respectively 
arc n + = \, ?i_=0; n + =l, n_=— 2; and n + =2, n_=l. 
So there are three resonant peaks appearing in optical 
absorption spectrum (see black curves in Fig. @] for 
cL=Q.21b), which correspond to the formation of mag- 
netoexciton states |2,1), |1,0), and |1,— 2), respectively, 
by absorbing photons of three different frequencies. As 
increasing the spacer thickness, these three absorption 
frequencies ux-^\f2(y/l-Q)=lAl, w 2 ^v / 2(l+V2)=3.41, 
and w 3 ^v / 2(\/2-l)=0.59 (see green dotted curves in 
Fig. 4). 

When the interlevel CI is also included, the above- 
discussed magnetoexciton states will be mixed and ener- 
getically redistributed, which observably correct the ex- 
act selection rule (ITU1) and bring about additional sub- 
sidiary peaks in the optical absorption spectrum. To 
clearly see this, we reconsider the case of /i=l as an ex- 
ample. The corresponding calculated optical absorption 
spectra for fully and partially occupied cases are respec- 
tively shown in Fig. [5ja) and[5jb). 

Let us first make an analysis on the fully occupied case. 
Because the interlevel CI mixes the noninteracting LLs, 
the absorption phenomenon should appear at the energy 
of every magnetoexciton state, with the resultant absorp- 
tion peak amplitude decided by the projection of the 
corresponding resonant excited states |n+,n_) allowed 



by the selection rule (T10|) in the case without interlevel 
interaction. When the spacer thickness d is small (e.g. 
d=0.2^s), not only the main resonant peak appears at 
the energy of excited level |2, 1), but also two weak res- 
onant peaks at the energies of excited levels |3, 1) and 
|4, 1) can be observed (see red curves in Fig. Ufa)). By 
increasing d to 2Mb, the resonant peak at |3, 1) become 




FIG. 5: (Color online) The optical absorption spectrums with 
[i=l for (a) fully and (b) partially occupied cases when the 
interlevel Coulomb interaction is included. The red, blue and 
green curves are for d=0.2ls, 2.01b, and 10.01b, respectively. 



in Fig. [Sfa)). When further increasing d to lO.OZs, we 
find that the resonant peak |3, 1) also disappears (green 
dotted curves in Fig. OJa)) and the optical absorption 
spectrum comes back to the noninteracting case. 

Similar physics also happens for partially occupied case 
(see Fig. E^b)). The overlap caused by the interlevel CI 
between magnetoexciton states will produce many sub- 
sidiary resonant peaks besides the main peaks at the en- 
ergies of magnetoexciton states |2,1), |1,0) and |1, — 2). 
Note that there is a subsidiary resonant peak appear- 
ing in front of |2, 1), which corresponds to the ground 
state |1,1). This special feature only appears for the 
partially occupied case and thus can be used to experi- 
mentally determine whether the highest occupied LL is 
fully filled or not in TIB. By increasing d, the subsidiary 
resonant peaks turn to disappear, and the optical ab- 
sorption spectrum turns back to the noninteracting case: 
only the main three resonant peaks at energies of states 
|2,1), 1 1,0) and |1,— 2) appear in absorption spectrum 
allowed by the selection rule (fTU|) . 

In summary, we have theoretically studied the opti- 
cal absorption properties of magnetoexcitons formed in 
TIB under a strong perpendicular magnetic field. By 
neglecting inter-Landau-level CI, we have determined an 
exact absorption selection rule, which greatly helps to 
connect absorption peaks with specific magnetoexcitons. 
The inclusion of inter-Landau-level CI has been shown to 
bring about subsidiary optical absorption peaks due to 
the mixing effect. Our present conclusion also works for 
graphene bilayers. 
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